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ABSTRACT 

We present a study of the X-ray to optical properties of a sample of 545 X-ray selected Type 1 AGN, from the XMM-COSMOS 
survey, over a wide range of redshifts (0.04 < z < 4.25) and X-ray luminosities (40.6 < Log Lp-iojkev ^ 45.3). About 60% of 
them are spectroscopically identified Type 1 AGN, while the others have a reliable photometric redshift and are classified as Type 1 
AGN on the basis of their multi-band Spectral Energy Distributions. We discuss the relationship between UV and X-ray luminosity, 
as parameterized by the a ax spectral slope, and its dependence on redshift and luminosity. We compare our findings with previous 
investigations of optically selected broad-line AGN (mostly from SDSS). A highly significant correlation between a ox and L 2500 a 
is found, in agreement with previous investigations of optically selected samples. We calculate bolometric corrections, £ bo i, for the 
whole sample using hard X-ray luminosities (Lp-iojkevX an d the Eddington ratios for a subsample of 150 objects for which black hole 
mass estimates are available. We confirm the trend of increasing bolometric correction with increasing Eddington ratio as proposed in 
previous works. A tight correlation is found between a ox and k\,„], which can be used to estimate accurate bolometric corrections using 
only optical and X-ray data. We find a significant correlation between a ox and Eddington ratio, in which a ox increases for increasing 
Eddington ratios. 

Key words, galaxies: active - galaxies: evolution - quasars: general - methods: statistical 



1. Introduction 

The distribution of X-ray to optical-UV ratio in quasars, as 
a function of optical and X-ray luminosity and redshift, has 
been the subject of activ e investigations for more than 30 years 
dTananbaum et al.f T979). The ratio is parameterized by the opti- 
cal to X-ray spectral index (a ox ) defined as: 



«ox = -- 



Log [L 2 kev/L 25 00A] 



2.605 



(1) 
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In the past, it was widely adopted to compute the QSO contri- 
bution to the X-ray background and esti mate the X-ray lumi- 
nosity function from optical counts (e. g., Zamo rani et alj fl981; 
lAvni & Tananbauiril Il986l ISilverman et alj 120051) . These stud- 
ies, based on Einstein observations, found a mean value of 
a ox in the range 1.3 — 1 .4 5. Subsequently, many investigations 
have been performed (e.g., [Yu an et al. 1998; lVignali e t al. 2003; 
Strateva et alj|2005t ISteffen et al|2006l hereafter S06 ; Uust et all 
20071 hereafter J0 7: iKellv et alj 120081: iGreen et all 120091 and 
Young et alll2009l) which were mainly based on AGN samples 
selected from the Sloan Digital Sky Survey (SDSS). Most of 
the previous studies on the relationship between a ox , luminos- 
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ity and redshift were based on large samples of optically se- 
lected type 1 AGN with a high fraction of X -ray detections. 
The S06 paper extends the lStrateva et al.l d2005l) work including 
moderate-luminosity optically selected AGN from the COMBO- 
17 survey with deep observations in the Chandra Deep Field- 
South. Their full sample consists of 333 optically selected Type 
1 AGN with a high fraction of X-ray detections (~ 88%) in 
the redshift range 0.01 < z < 6.28. More recently, J07 have 
further extended the S06 sample with 59 optically luminous 
quasars spanning a red shift range of 1 .5-4.6 and 14 quasars from 
IShemmer et all j|2006l) . All the X -ray information on their sam- 
ple is derived from Chandra, XMM-Newton and ROSAT ob- 
servations. The X-ray to optical flux ratio distribution of an X- 
ray selected sample of 188 Type 1 AGN from the ChaMP sur- 
vey, which is limited to bright X-ray fluxes and optical lumi- 
nosities (f x > 10 14 erg cm^s' 1 and r < 22.5), is presented in 
ISilverman et alj d2005l) . A strong correlation between a ox and 
the optical luminosity at 2500A is foun d, while a ox is only 
marginally dependent upon redshift (but see Bec htold et al.l2003l 
for different results). The a ox distributions typically cover the 
range 1.2-1.8, with a mean value of about 1.5. 

These studies could not properly address the effect of a 
source selection in a different band (i.e. X-rays). The availabil- 
ity of large samples of X-ray selected QSOs with a high quality 
photometric and spectroscopic coverage in the optical now opens 
the possibility for an extended investigation of the a ox distribu- 
tion and its evolution. 

Understanding how a ox evolves with luminosity and redshift 
may provide a first hint about the nature of the energy generation 
mechanism in AGN. It is also a first step to wards an estimate of 
the AGN bolometric luminosity function (Hop kins et alj|2007h 
and t he mass function of Supermassive Black Holes (SMBHs) 
(e.g. lMarconi et all f2004). and towards the understanding of the 
structure of the AGN accretion disk and X-ray corona. While the 
calculation of the bolometric luminosity requires good-quality 
data over a large portion of the electromagnetic spectrum, it is 
relatively easy to compute a ox for sizable samples of objects up 
to high redshift. 

We analyze the dependence of a ox upon redshift, optical 
and X-ray luminosities using a large X-ray selected sample 
of Ty pe 1 AGN in t he Co smic Evolution Survey (COSMOS) 
field dScoviHe et alj 120071) . The COSMOS field is a so far 
unique area for deep and wide comprehensive multi-wavelength 
coverage: radio with the VLA, infrared with Spitzer, optical 
bands with Hubble, Subaru, SDSS and other ground-based 
telescopes, near- and far-ultraviolet bands with the Galaxy 
Evolution Explorer (GALEX) and X-rays with XMM-Newton 
and Chandra. The spectroscopic coverage with VIMOS/VLT 
and IMACS/Magellan, coupled with the reliable photometric 
redshifts derived from multiband fitting, allows us to build a 
large and homogeneous sample of QSOs with a well sampled 
spectral coverage and to keep selection effects under control. 

The broad-band information contained in the COSMOS 
database is well suited for a detailed study of AGN Spectral 
Energy Distributions (SEDs), bolometric luminosities (Lboi) and 
bolometric corrections, in particular the one from the X-rays, 
defined as: 

^bol = Lbol/Lp-KTjkeV- (2) 

If the black hole mass is available (Mbh), it is also possi- 
ble to compute the Eddington ratio (/lEdd = Lboi/LEdd, where 
LEdd = 1-38 x 1O 38 Mbh/M ) and investigate any possible de- 
pendence of /lEdd on the SED shape. In particular, it is possible 
to test whether the correlation between hard X-ray bolometric 



Table 1. Selection criteria from 1312-XMM catalog. 



Catalog 


BL" 


Photoz 


Total b 


1312XMM C 


361 


613* 


974 


Class > 19 d 


(361) f 


236 


597 


Radio quiet c 


322 


223 


545 



a Broad-line sources from optical spectroscopy (FWHM >2000 km 

b This column gives the sum of the objects with spectroscopic or pho- 
tometric redshifts. 

c F[o.5-2]kev - 10~ 15 erg sec^'cnT 2 and secure optical association. 
The 1312-XMM catalog also comprises: 241 sources which are not 
broad-line AGN; 34 sources which are spectroscopically classified 
as stars; 49 sources which are best-fitted with a stellar template and 
14 sources for which th e photometric redsh ift is not reliable. 

d According to Table 2 in lSalvato etafl i2009h . 

e if < 10. Radioloudness is defined in eq. {3). We exclude 51 radio- 
loud sources: 39 with spectroscopic redshift and 12 with photomet- 
ric redshift. 

f All spectroscopically identified broad-line AGN are included, even 

if the classification is less than 19. 
* Only extragalactic sources are considered. 



correction fcboi and Eddington ratio /lEdd, recently obtained for 
bright AGN in the nearby Universe ( Vasudevan & Fabian 2009; 
hereafter VF09), holds also at higher redshifts. 

This paper is organized as follows. In Section[2]we describe 
the selection criteria for the sample used in this work. Section [3] 
presents the data and the method by which we construct SEDs. 
Data analysis and results are given in Section[4] In Section|5]we 
estimate the possible effects of reddening and host-galaxy light 
contribution on our main results. The discussion of our findings 
is given in Section [6] and a summary of the results is given in 
Section|7] We adopted a flat model of the universe with a Hubble 
constant Hp - 70 km sec" 1 Mpc" 1 , Q.m = 0.29, — 1 - D.m 
dKomatsu et al.ll2009t) . 

2. The Data Set 

2. 1 . The Parent Sample 

The XMM-COSMOS catalog comprises 1822 point-like X-ray 
sources detected by XMM-Newton over an area of ~ 2deg 2 
for a total of ~ 1.5 Ms at a fairly homogeneou s depth of 
~ 50 ks dHasinger et al.ll2007llCappeUuti et al.ll2009l) . Following 
iBrusa et al.l d2009b ), we excluded 24 sources which turned out 



to be a blend of two Chandra sources leading to a total of 1798 
X-ray selected point-like sources. 

Spectrosc opic redshif t s for th e proposed counterparts are 
compiled by IBrusa et all d2009bl) from the Mag ellan/IMACS 
and M MT observation campaigns (~ 590 objects, Trum p et al 
2009;), from the zCOSMOS project (~ 350 objects. iLillv et al 
120071) . or were already p resent either in the SDSS survey 
catalog (~ 100 objec ts, lAdelman-McCarthv & et all 120051 
i Kauffmann etail 120 03?!). or in the literature (~ 95 objects, 



Prescott et al.l |2006). In summary, good-quality spectroscopic 



redshifts are available for 738 sources, corresponding to a sub- 
stantial fraction (~ 45%) of the entire XMM-Newton sample. 



1 These sources have been retrieved from the Nasa Extragalactic 
Database (NED) and from the SDSS archive. 
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Photometric redshifts for almost all XMM-COSMOS 
sources have been obtained exploiting the COSMOS multi- 
wavelength database and are presented in Salvato et al. (2009, 
hereafter S09). Since the large majority of the XMM-COSMOS 
sources are AGN, in addition to the standard photometric red- 
shift treatments for normal galaxies, a new set of SED tem- 
plates has been adopted, together with a correction for long- 
term variability and luminosity priors for point-like sources (see 
below and S09 for further det ails). The availability of the in- 
termediate band Subaru filters (Tanig uchi et al.ll2007l) is c rucial 



Galaxy-Type 2 SED 



Type 1 SED 



in picking up emission lines (see also Wolf et al. 2004). This 



led, for the first time for an AGN sample, to a photometric red- 
shift accuracy comparable to that achieved for inactive galax- 
ies (cta z / ( i + ,) ~ 0.015 and ~ 4% outliers) down to i ^22.5. At 
fainter magnitudes (22.5 < i < 24.5), the dispersion increases to 
°"Az/(i+z) - 0.035 with ~ 4.8% outliers, still remarkably good for 
an AGN sample. A photometric redshift is available for all but 
32 objects out of 1798. 

In addition to the photometric redshifts, S09 provide also a 
photometric classification based on the best-fit broad-band SED 
template. Briefly, each AGN SED has been fitted with a total 
of 30 different templates which include both normal galaxies 
(early type, late type and ULIRG galaxies), low-and-high lumi- 
nosity QSOs (both Type 1 and 2) and hybrids created assum- 
ing a varying ratio between the AGN and a gal axy templates 
(90:1 0, 80:20,..., 10:90; see S09 for details and iPolletta et all 
120071) . About 40% of the sources are best-fitted by AGN- 
dominated SED, while the remaining sources are reproduced by 
host galaxy-dominated SED. The photometric classification is 
also confirmed a posteriori (see Fig. 10 in S09) with the distri- 
bution of the XMM sources in the X-ray h ardness ratio plane 
dCappelluti et al.ll2009HHasinger et alj|2007b . 

We restricted the analysis to the X-ray sources detected 
in the soft band at a flux larger than 10~ 15 erg sec _1 cm~ 2 
(Cappelluti et al. 2009). Given that objects for which no secure 
optical counterpart could be assigned, are usually affected by 
severe blending problems, making the photo-z estimate unre- 
liable, our parent sample consists of 1312 sources (hereafter 
1312-XMM) for which a secure optical counte rpart can be as- 
sociated (see discussion in lBrusa et al.ll2009blal) . 



2.2. Type 1 AGN Sample 

From the 1312-XMM catalog we have selected 361 spectroscop- 
ically classified broad-line AGN^on the basis of broad emission 
lines (FWHM > 2000 km sec -1 ) in their optical spectra. We will 
refer to this sample as the "spectro-z" sample. As a comparison, 
in the 1312-XMM catalog there are 241 objects spectroscopi- 
cally classified as not broad-line AGN (Type 2 or emission-line 
or absorption-line galaxies). 

The spectroscopic completeness of the X-ray selected sam- 
ple quickly decreases towards faint optical magnitudes. A siz- 
able fraction of Type 1 AGN may be present among optically 
faint X-ray sources; not to use them in our analysis would sig- 
nificantly bias the distribution of X-ray to optical flux ratios. In 
order to extend our Type 1 AGN sample to fainter magnitude, 
we proceed as follows. 



2 The origin of spectroscopic redshifts for the 361 sources is as 
follo ws: 63 objects from the SDSS archive, 75 from MMT observa- 
tions dPrescott et alj2006h . 1 12 from the IMACS observation campaign 
(lTrumpet alJl2007h . 93 from the zCOSMOS bright 10k sample (see 
iLillv et al.12 007) and 18 from the zCOSMOS faint catalog. 
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Fig. 1 . Distribution o f the photometric classification (see Table 2 
in lSalvato et al.ll2009b for the spectroscopically identified Type 
1 AGN sample (hatched histogram), for the photometric sample 
(open histogram) and for the not-Type 1 AGN sample (filled his- 
togram). The dashed line at SED Type=19 marks our adopted 
separation between sources with photometric redshift that we 
define as Type 1 AGN (236 sources on the right) and not-Type 
1 (377 sources on the left) on the basis of the best-fitting SED 
template. 



• 361 spectro-z. 
x 236 photo-z 




Fig. 2. Plot of the [0.5 - 2]keV flux versus the total f CFHT 
magnitude for the spectroscopically identified (black points) 
and the photo-z classified (red crosses) broad line AGN. The 
dashed lines represent the constant X-ray to optical flux ratio 
log(X/0) = ±1. 
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Fig. 3. Redshift distribution of the 545 Type 1 AGN consid- 
ered in this work (open histogram). The filled histogram shows 
the redshift distribution for the sample of 322 spectroscopically 
identified sources, while the hatched histogram is the redshift 
distribution for the 223 sources without spectroscopic redshift. 

First, we looked at the distribution of the best-fit photometric 
classifications for the spectroscopically identified sample. Figure 
Q] shows the distribution separately for the 361 spectroscopi- 
cally identified Type 1 AGN (hatched histogram) and for the 
241 sources which are not broad-line AGN (filled histogram). 
The dashed line at the SED Type=19 marks the division be- 
tween galaxy-dominated and AGN-dominated SED (see Table 2 
in S09 for details). The large majority of the broad emission line 
AGN in the spectro-z sample (> 90%) are classified as Type 1 
AGN by the SED fitting. The distribution of the 241 not broad- 
line sources is in agreement with the SED-based classification 
(~ 83% have Class<19), and the number of not-Type 1 which 
have Class>19 is relatively small (42 sources, about 17%). Then, 
we consider the remaining 613 X-ray sources in the 1312-XMM 
sample for which only a photometric redshift is available from 
S09, and exclude all the sources which are best-fitted with a stel- 
lar template. The distribution of the best-fit templates for this 
photometric sample is also shown in Fig. Q] (open histogram). 
In the following, we assume that the 236 X-ray sources, classi- 
fied by the SED fitting with an AGN-dominated SED are Type 1 
AGN. We will refer to this sample as the "photo-z" sample. 

In Figure [2] we plot the soft X-ray flux as a function of i* 
CFHT magnitude. Points and crosses represent sources in the 
spectro-z and photo-z sample, respectively. The dashed lines de- 
limit the region typically occupied by AGN along the X-ray to 
optical flux ratio \og(X/ O) — +10 The 236 sources in the photo- 
z sample increase the completeness of the data-set at fainter 
magnitudes. While some Type 1 AGN may have been missed 
(about 10%, with SED Type<19 in Fig. Q]), we are confident 
that the described choice minimizes the selection bias against 
optically faint Type 1 AGN. We only include sources that have 
optical-UV data in order to compute monochromatic luminosi- 

3 \og(X/0) = log/,. + f/2.5 + 5.6 



ties at 2500A and the SED. This selection criterium only ex- 
cludes the photometric source with XID=5120. 

It is well known that radio-loud AGN have an enhanced X- 
ray emission mechanism linked to the jets, which can provide an 
increment in the X-ray emission with resp ect to radio-quiet AGN 
with similar optical luminosities (e.g. , IZamorani et alj 1 198 lL 
IWilkes & Elvisfl 19871 ICappi et al.lfl997h . We exclude RL AGN 
from the total sam ple using the "standard" definition of radio- 
loudness, R > 10 rtKellerma nnet alJfl989 ; )Pl We converted the 
monochromatic flux at 1 .4 GHz, reported in the final catalog of 
the VLA-COSMOS Deep project (see for mo re detail Schinnerer 
et al. 2009, submitted, and B ondi et alj20 08) to L5ghz(v) assum- 
ing f(v) oc v~ a with a = 0.7. We excluded from the sample 51 
radio-loud Type 1 AGN (39 AGN from the spectroscopic sam- 
ple and 12 sources with photometric redshift) with a value of 
radio-loudness R > 10. 

The final Type 1 AGN sample used in our analysis, there- 
fore, comprises 545 X-ray selected AGN (322 from the spectro- 
z sample and 223 from the photo-z sample) spanning a wide 
range of redshifts (0.04 < z < 4.25) and X-ray luminosities 
(40.6 < Log L[2-iojkeV ^ 45.3). The selection criteria are sum- 
marized in Table Q] 

Assuming that the fraction of misclassificated sources in 
the spectro-z sample (17%) and the fraction of missing Type 1 
sources (10%) could be applied to the photometric sample, we 
are able to estimate the uncertainties associated to the photo-z 
sample. If this were the case, the incompleteness and the con- 
tamination on the total sample are very low, namely about 4% 
and 7%, respectively. 

The redshift distributions of the total, spectroscopic and pho- 
tometric samples are presented in Figure[3] The median redshift 
of the total sample is 1.57 (the mean redshift is 1.66, with a dis- 
persion of 0.75). The median redshift of the spectro-z sample is 
1 .54, while the median of the photo-z sample is 1 .66. 

3. Rest-frame monochromatic fluxes and 
Spectral Energy Distributions 

To obtain rest-frame monochromatic luminosities at 2 keV and 
2500A and estimate bolometric luminosities, we used all the 
multi-color information as compiled by Brusa et al (2009). The 
catalog includes multi-wavelength data from mid infrared to 
hard X-rays: MIPS 24 pm G03 da ta (Le Floc'h et al 2009, 
ApJ submitted), IRAC flux densities (Sanders et al. 2007), near- 
infrared CFHT/K-band data dMcCracken et al.l2008l). HS T/ACS 
F814W imaging of the COSMOS field dKoekemoer et al.ll2007l) 
optica l multiband photometry (SPSS. Subaru. ICapak et all 
120071) and near- and far-ultraviolet bands with GALEX 
(Zamojski et al. 2007). More specifically, the number of detec- 
tions at 24 pm is 472; for the 73 undetected sources, we consider 
5cr upper limits of 0.08 mJy. Most of the sources are detected 
by IRAC: 545, 543, 544 and 543 at 3.6 pm, 4.5pm, 5.8 pm and 
8.0pm band, respectively (all Type 1 AGN were detected in the 
3.6pm IRAC band). For the undetected sources we consider 5<x 
upper limits of 1.7 pjy, 1 1.3 p}y and 14.6 p Jy at 4.5 pm, 5.8 pm 
and 8.0pm, respectively (see Table 3 of ISanders et al.l 12007). 

4 Radio-loudness is defined as 



where L$ghz(v) and Lb(v) are the rest-frame monochromatic luminosi- 
ties at 5 GHz and in the optical B band, respectively. 
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Fig. 4. The SED of a spectroscopically identified QSO atz = 1.8 
(XID=2046, black line), compared with the mean SED of Elvis 
et al. (1994, blue line). The rest-frame data, used to construct 
the SED of XID=2046 are, from left to right: black triangle, 
24 fim MIPS; red square, 8.0ji/m IRAC; green square, 5.8 yt/m 
IRAC; cyan square, 4.5 ymn IRAC; blue square, 3.6 yum IRAC; 
grey hexagon, K CFHT; blue triangle, J UKIRT; black, red, 
green, cyan and blue hexagons represent z + , r + , g + , Vj and Bj 
SUBARU bands, respectively; black and orange squares repre- 
sent i* and u* CFHT bands, respectively; black hexagons soft 
X-ray and medium-hard X-ray luminosity. 



Only very small objects went undetected in the optical and near 
infrared bands: only 2 upper limits in the z + band; 1 upper limit 
in the B j, Vj and r + bands; 8 upper limits in both i* and u* CFHT 
bands; 9 in the K s CFHT band and 31 in the J UKIRT band. 
The observations are not simultaneous, as they span a time in- 
terval of about 5 years: 2001 (SDSS), 2004 (Subaru and CFHT) 
and 2006 (IRAC). In order to reduce variability effect, we se- 
lected the bands closest in time to the IRAC observations (i.e., 
we excluded SDSS data). In Table|2]we list the main X-ray and 
optical properties of the sample. The data for the SED computa- 
tion were blueshifted to the rest-frame and no K-correction has 
been applied. Galactic reddening has been taken into account: 
we used the selective at tenuation of the st ellar continuum k(A) 
taken from Table 11 of Capaketal. ( 20071) . Galactic extinction 
is estimated from ISchlegel et alj d!998l) for each object in the 
1312-XMM catalog. Count rates in the 0.5-2 keV and 2-10 keV 
are converted into monochromatic X-ray fluxes in the observed 
frame at 1 and 4 keV, respective ly, using a Galactic co lumn den- 
sity N H = 2.5 x 10 2l) cm~ 2 (see lCappelluti et alj|2009l) , and as- 
suming a photon index F x = 2 and F t = 1 .7, for the soft and hard 
band, respectively. 



3.1. The SED Computation 

With a procedure similar to that adopted by Elvis et al. (2009, 
in preparation), we assembled SEDs for the 545 Type 1 AGN. 
First of all, we plotted the information of all the sources, from 



Fig. 5. The SED of a spectroscopically identified QSO at z 
1.161 (XID=265, black line). Keys are as in Fig. |4] 



24 yum to the hard X-ray band (the few upper limits in the optical 
and near-infrared are not considered), in the Log v - Log [vL v ] 
rest-frame plane. For each source we consider the rest-frame lu- 
minosity and corresponding frequency in each of the available 
bands. We determine a first order polynomial between two adja- 
cent points, and the resulting function is then sampled with a step 
of ALog v = 0.085. In this way a "first order" SED, where ev- 
ery point is separated by 0.085 in log frequency, is created. The 
resulting SED is "smoothed" using a fourth order polynomial in- 
terpolation and re-sampled with a step of ALog v = 0.01. This 
choice is motivated by the fact that a single interpolation with 
a high-order polynomium could introduce spurious features in 
the final SED. This procedure allows us to build densely sam- 
pled SEDs at all frequencies, to extract accurate 2500 A and 2 
keV rest-frame monochromatic luminosities at all redshifts, to 
estimate bolometric luminosities and bolometric corrections. 

For the computation of the bolometric luminosity we need 
to extrapolate the data in the UV to X-ray gap and at high X- 
ray energies. From the rest-frame UV luminosity data point at 
the highest frequency in each SED, we assume a power law 
spectrum to 500A, as meas ured by HST obse rvations for radio- 
quiet AGN (f v oc v -1 - 8 , see lZhenget al.lll997l) . We then linearly 
connect the UV luminosity at 500 A to the luminosity corre- 
sponding to the frequency of 1 keV. We extrapolate from the 
X-ray data point to the luminosity at 1 keV computing the slope, 
a x = Log (f4kev/fikev)/Log (v 4keV /vi keV ) (a x = T x - 1). Finally, 
the X-ray spectrum is extrapolated at higher energies using the 
observed slope a x , and introd ucing an expon ential cut-off at 200 
keV, f(E) ~ £-«v e -£/200keV jGim et alj|2007b . Hence, the bolo- 
metric luminosity is computed integrating the extrapolated SED 
in the Log v - Log [vL v ] rest-frame plane. Our bolometric lu- 
minosities are in agreement with those computed by Elvis et al. 
(2009). In Figure|4]and|5]we show two examples of SED used for 
the calculation of the bolometric luminosities and the rest-frame 
monochromatic luminosities of the Type 1 AGN sample. Given 
the selection criteria, the SEDs of the objects in our sample are 
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Table 3. Correlations and their significance from Kendall-r 
statistics (K-t) and from Partial Kendall-r. 



Correlation 


K-T 


Controlling variable 


Partial K-t 




(0-) 




(cr) 


Total 


Q'ox _ L 2500 ^ 


17 


z 


21 


Q'ox _ L 2 keV 


1.2 


z 


1.5 


a: ok - Z 


7 


L2500 A 


< 0.1 


Spectro-z 


«ox - L 25 oo A 


12 


z 


14 


ttox - L2 keV 


0.2 


z 


2.3 


•^ox - Z 


5 


L2500 A 


< 0.1 


Photo-z 


«ox - L 25 oo A 


10 


z 


12 


Q'ox _ L 2 keV 


0.2 


z 


3 


«ox - Z 


5 


L2500 A 


0.7 



quite differe nt, ranging from b lue QSO with an SED close to that 
reported by lElvis et all (1 19941) for radio-quiet AGN, to objects 
with red SED possibly due to intrinsic absorption and/or host- 
galaxy contamination. We will discuss in Sections 15.11 and 15.21 
how we take the possible effects of reddening and host-galaxy 
contamination into account in our analysis. A complete Type 1 
AGN SED atlas will be presented in Elvis et al. (2009). 



4. Statistical Analysis 

To study the possible correlations between a ox , L 250oA , L2keV 
and redshift, we used the Ast ronomy Survival Analysis soft- 
ware package (ASURV rev. 1.2; IIsobe et alJl990tlLavallev et alj 
Il992l) . ASURV implements the bivariate data-analysis methods 
and also p roperly treats censored data using the surviv al analysis 
methods dFeigelson & NeIsonlll985l; llsobe et alJ[T98 6). We have 
employed both full parametric estimate and maximized (EM) re- 
gression algorithm and semiparametric Buckley-James regres- 
sion algorithm (Buckley & James 1979) to perform the linear 
regression of the data. The EM regression algorithm is based on 
the ordinary least-squares regression of the dependent variable 
Y against the independent variable X (OLS[Y|X]). The regres- 
sion line is defined in such a way that it minimizes the sum of the 
squares of the Y residuals. Since in our analysis the choice of the 
independent variable is not straightforward, we have also used 
the inverse of OLS(Y|X) (OLS[X|Y]) and we then calculate d 
the bisector of the two regression lines (see llsobe et al.| [i990). 
We only present the findings from the EM regression algorithm, 
since in all cases the results from the Buckley-James regression 
algorithm agreed within the errors; we also report the findings 
from OLS bisector. Moreover, we used the partial-correlation 
analysis method to compute the correlation between two vari- 
ables, checking the effect of an additional parameter which the 
two variables depend on. This method is implemented in the 
FORTRAN program PARTIAL.TAU, available from the Penn 
State Ce nter of Astrostatist i cs, an d uses the methodology pre- 
sented in Akritas & Sieberj dl996l) . The procedure is based on 
Kendall r-statistic that properly handles censored data. 

We summarize the results for partial-correlation analysis in 
Tabled 



41 ■ L 250()A VsL 2keV 

Previous studies on optically selected AGN reported 
a relationship between L2k e v and L 250()A in the form 
L2kev x ^2500 A' w ' tn best-fit values of the exponent 
between 0.7 H - 0.8 (e.g. lAvni & Tananbauri fl98l 1986: 
Chanan 1983: iKriss & Canizare sj 119851: lAnder son & Margon 
1981 IWilkes et al.lll994t lYuan et al.lll998UVignali et al.l l2003b 
Stratevaetal] 120051: S06; J07; see also lLa Franca et all d 1995b . 
who found f3 ~ 1). For the best-fit parameters using OLS(Y|X) 
(i.e. treating L 2500 A as the independent variable) we find 

Log L 2 keV = (0.599 + 0.027)Log L 2500 A + (8.275 ± 0.801) (4) 

while treating L 2500 A as the dependent variable (i.e., OLS(X| Y)) 
and inverting the resulting best-fit, we find 

LogL 2k ev = (0.952 + 0.033)LogL 250oA -(2.138+0.975). (5) 

We then c ompute the bisecto r of the two regression lines as de- 
scribed by llsobe et alj(fl9 90) and find 

Log L 2 keV = (0.760 + 0.022)Log L 2500 A + (3.508 ± 0.641) (6) 

with a dispersion of 0.37. The difference of the best-fit f3 with 
respect to a linear correlation (i.e., f3 = 1) is highly signifi- 
cant (11 cr), this result confirms the non linear correlation be- 
tween L 2 kev - L 250oA . This implies not only that a ox must be 
dependent on optical luminosity, but also that optically bright 
AGN emit less X-ray (per unit UV luminosity) than optically 
faint AGN. Comparing our results with those obtained from opti- 
cally selected samples, we find that our slope is consistent within 
2 cr with those of S06 sample (J3 S06 = 0.721 + 0.011) and J07 
sample (J3jcn = 0.709 + 0.010), while it is significantly smalle r 
than the slope of 1.117 + 0.017 found by iGreen et all (12 009). 
Moreo ver, our slope is similar to the value found bv lStalin et all 
(120091) who used an X-ray selected sample composed by 267 
broad-line AGN (J3 = 0.870). Treating separately the spectro-z 
and photo-z samples, we found that the best-fit slope for the bi- 
sector is 0.782 ± 0.033 with a normalization 2.815 + 0.989 and 
0.786 + 0.033 with a normalization 2.824 + 0.958, respectively. 
The spectro-z sample is limited to relatively bright 2500 A lumi- 
nosities, while the photo-z sample extends the luminosity range 
of about 1 .5 dex towards lower L 2500 A values. It is worth noting 
that there is a significant difference between the normalization 
in eq. (O and the normalizations computed treating separately 
the spectro-z and the photo-z samples. The percentual difference 
between the predicted L 2 k e v using the best-fit relation for the 
total sample and the best-fit for the spectro-z sample is ~ 10%. 
Moreover, considering only the results from the spectro-z sam- 
ple, we would slightly underestimate the dispersion in the bi- 
sector. These considerations justified the choice of adding the 
photo-z sample in order to reduce possible biases associated to 
the spectro-z sample. In Figure [6] we show monochromatic X- 
ray luminosities as a function of monochromatic UV luminosi- 
ties for our X-ray selected sample. For comparison, we plot the 
L2 kev - L 2500 A relations from S06 and J07. In the bottom panel, 
we report residuals from equation while on the right panel, 
residuals from eq. © with the mean and the lcr dispersion com- 
puted for each bin with ALog L = 1 are reported. 

4.2. a ox i/sL 2500/ 4 

Previ ous X-ray studies of AGN using optical ly selected samples 
(e.g.: IVignali et al.ll2003LlStrateva et al.l2005l S06, J07) revealed 
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Table 2. Optical and X-ray properties of the total sample. 



XID 


Spectroscopic 


Photometric 


kog ^2500 A 


Log L 2 kcv 




Log L[2-10]keV 


Lbol,l^m 


kbol 


Class" 




redshift 


redshift 


[erg s^'Hz -1 ] 


[erg s~'Hz~'] 




[ergs -1 ] 








1 


0.373 


37 +0 - 03 
-fl.fll 


29.49 


26.04 


1.32 


44.00 


11.53 


12.94 


28 


2 


1.024 


1 s-* r- _i_n ni 

1 05 


29.86 


26.95 


1.11 


44.97 


12.30 


8.09 


30 


5 


0.345 


n 1 £+0 02 

0.36 


29.22 


26. UO 


1 1A 

1.24 


a i no 
43.98 


1 1.44 


11 K 


29 


6 


0.360 


32 +m 

"• J -0.02 


29.22 


25.43 


1.45 


43.39 


11.32 


32.60 


20 


8 


0.699 


071 +o.oi 
0.68+jl 


30.00 


26.27 


1.43 


44.09 


11.98 


29.60 


26 


10 


0.689 


29.24 


25.92 


1.28 


43.71 


11.30 


14.93 


22 


13 


0.850 


0.84!| 


29.54 


26.06 


1.34 


43.98 


11.70 


20.06 


24 


15 


2.033 


2 - 01 -l 
061+ 8.«§ 

-QQl 


30.83 


27.11 


1.43 


45.11 


12.95 


26.33 


30 


16 


0.667 


29.32 


25.82 


1.34 


43.78 


11.38 


15.14 


19 


17 


1.236 


1 oo +0.02 


30.77 


26.77 


1.53 


44.70 


12.75 


43.09 


29 



Notes — This table is presented entirely in the electronic edition; a portion is shown here for guidance. 



a The SED model of the best-fit template is coded from 1 to 30, as detailed in S09. 



2 - 



_This work 
. S06 




28 30 

Iog[L(l>) a500 A ] [erg s-> Hz"'] 



32 



Fig. 7. Plot of the a ox vs. the rest-frame monochromatic lu- 
minosity L 



2500 A 



for our X-ray selected sample: 322 spectro- 
scopic sources (circles) and 223 sources with photometric red- 
shift (crosses). The solid line represents our best-fit given by 
equation (0. For comparison, the best-fit derived by S06 (short- 
dashed line) and from J07 (long-dashed line) are shown. The 
lower panel shows the residuals (Aa ox ) of our best-fit relation. 
The error bars represent the mean and the lcr standard deviation 
of the mean of the residuals for each ALog (L 2500 \) = 1 bin. 



a highly significant correlation between a ox and the 2500 A 
monochromatic luminosity, with a slope of ~ 0.14. We apply 
EM regression to our X-ray selected Type 1 sample and confirm 
the a ox - L 2500 \ correlation at the 17 cr significance level. The 
correlation becomes stronger when taking into account the effect 
of redshift using Kendail-r partial-correlation analysis (~ 21 cr). 
The best-fit relation for a ox - L 2500 \, using OLS(Y|X) (i.e. treat- 
ing L 2500 A as the independent variable), is 

«ox(L 250oA ) = (0.154+0.010)LogL 250oA -(3. 176+0.223), (7) 



with a dispersion of 0. 18. For the spectro-z sample, we found a 
slope of 0.166 + 0.011 with a normalization of -3.541 + 0.335, 
while for the photo-z sample the best-fit slope is 0.142 + 0.012 
with a normalization -2.83 1 +0.345. In Figure|7]we plot a ox ver- 
sus L 250oA for our sample. For comparison, we plot the best-fit 
linear regression found in S06 (short-dashed line) and J07 (long- 
dashed line). Our slope is statistically consistent within ~ 1.5<x 
with the slopes of the S06 and J07 best-fits. Moreover, our 
slope is c onsistent, within 1 .6c and lcr, with those published re- 
cently bv lKellv et ail d2008h (0. 1 2+0.02) andl Young et all d2009h 
(0. 1 53 + 0.0 1 2), res pectively. Differently. iGreen et alJ d2009l) and 
IStalin et al.l d2009t) found a flatter slope of 0.061 ± 0.009 and 
0.065 ± 0.019, respectively. The residual Aa ox , defined as 



ffoxCL^An a) 



(8) 



where a ox values are measured from the observations, are shown 
in the bottom panel of Fig. [7] The bars represent the mean and 
the 1 cr dispersion for ALog (L 2500 A ) = 1 bins. 

4.3. or x vs L2 kev and vs redshift 

We find no significant correlation between q ox and L2k e y (e-g-i 
iRriss & C anizares 1985: Avni & Tanan baum 1986t lWilkes et al.l 
1994 and Yuan et aLlll998l) . For comparison with S06 and J07, 
we perform the EM regression method and the best-fit parame- 
ters for the a ox - L2k e v relation (i.e. treating L2kev as the inde- 
pendent variable) are 

«ox(L 2 kev) = (0.019 + 0.013)Log L 2 ke v + (0.863 + 0.344). (9) 

For the spectro-z sample the slope is 0.012+0.020, while for the 
photo-z sample it is -0.020 + 0.019. Analyzing in the same way 
the a ox - z relation (i.e. treating z as the independent variable), 
we find 



a ox ( z ) = (0.064 + 0.010) z + (1.261 + 0.018). 



(10) 



The significance of the correlation between a ox and redshift is 
7 cr. This may suggest a not negligible evolution of a ox with red- 
shift, but if we account for the effect of the optical luminosity 
using partial Kendall-r statistics, the significance of the correla- 
tion disappears (< 0.1 cr). In Figure [8] and [9] we show a ox as a 
function of L2 kev and redshift, respectively, where the solid line 
represents the best-fit relation (eq. J9) and ifTUl ). In the bottom 
panel of each plot we show the residuals and the lcr dispersion 
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log[L(i/) 2500 A ] [erg s" 1 Hz" 1 ] 

Fig. 6. Values of the rest-frame monochromatic luminosity L2 kev plotted versus the rest-frame L 2500 - K monochromatic luminosity for 
our X-ray selected sample: 322 spectroscopic sources {circles) and 223 sources with photometric redshift (crosses). The solid black 
line shows the best-fit relation that we found using the OLS bisector algorithm (see eq. 0). For comparison, the best-fit derived 
by S06 (short-dashed line) and the best-fit from J07 (long-dashed line) are also shown. The orange and the green line represent 
equations (O and Q, respectively. The lower panel and the panel on the right show the residuals of log L2 kev and log L 2500 ^ with 
respect to the best-fit lines given by eq. H) and |5)), respectively. The error bars represent the mean and the Icr dispersion of the 
residuals for each ALog (L 25()0 s ) = 1 bin. 



in bins of Az = ALog L2 kev = 1 ■ The absence of a correlation 
between a ox and L2 kev is consistent with the fact that the cor- 
relation between L 25on x and L2kev, treating L2k e v as the inde- 
pendent variable (eq.jS)), is very close to be linear (i.e., /3 ~ 1). 
From the partial correlation analysis we showed that there is no 
correlation between a ox and redshift, once the effect of optical 
luminosity is properly removed. The energy mechanisms that 
generate the broad-band emission in AGN do not vary, in an ap- 
preciable way, over cosmi c time. This is in a gree ment with most 
of prev ious works (but see I Yuan et al.l dl998) and lBechtold et al.l 
(2003) who claimed that the primary correlation comes out upon 
redshift). The difference in the normalization between our ff ox -z 
best-fit relation and S06 best-fit relation is probably due to the 
different optical luminosity range spanned by our data and S06 
data. The average optical luminosity of the S06 sample is more 
than 1 decade higher than ours. Because of the correlation be- 
tween a ox and L 2 5 oA> m i s implies a stepeer a ox at any given 
redshift. For the spectro-z sample we found that the slope is 
0.058 + 0.012, while for the photo-z sample it is 0.079 + 0.016. 



4.4. Bolometric Corrections vs Eddington ratio 

In Figure [TO] we plot a ox vs. k\, Q \. In order to better estimate the 
bolometric correction values we select a subsample of 343 AGN 
with detections in both soft and hard bands. The best-fit relation 
is computed using a second-degree polynomial: 

Log k bol = 1.561 - 1.853a ox + 1.226o4- (H) 

We quantify the lcr (~ 0.078) and the 3<x (~ 0.234) dispersion in 
the Log kboi-ffox relation using a 3.5<x clipping method. A linear 
correlation is not a good description of the observed data points. 
The addition of a quadratic term significantly improves the fit 
quality (A% 2 ~ 140, significant at > 8<x confidence level ac- 
cording to an F-test). A visual inspection of the individual SED 
of the sources which mostly deviate from a linear fit (i.e. those 
with a ox < 1.2) indicate that their SEDs are not significantly 
different from those of the rest of the sample. We also note that 
the dispersion around the best-fit is very small For 150 Type 1 
AGN in our sample in the redshift range 0.196 < z < 4.251, 
an estimate of the BH mass is available from virial estimators 
(Peterson et al. 2004) using the Mg n line width (63 sources from 
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Fig. 8. Plot of a ox vs. the rest-frame monochromatic luminosity 
L2kev for our X-ray selected sample: 322 spectroscopic sources 
(circles) and 223 sources with photometric redshift (crosses). 
The solid line represents our best-fit given by equation (|9). For 
comparison, the best-fit derived by S06 (short-dashed line) and 
from J07 (long-dashed line) are also shown. The lower panel 
shows residuals (Aff ox ) from our best-fit relation, while error 
bars represent the mean and the lcr dispersion of the residuals 
for each A[Log L,2kev] = 1 bin. 



Merloni et al. submitted; 63 sour ces from Trump et al.l2009f) : the 
Hj6 line width (16 sources fromlTrump et alj|2009h and the C iv 
line width (8 sources from Trump et al. 2009|). 

The bolometric luminosity which enters in the calculation of 
the Eddington ratio, /l E dd, is computed by integrating the individ- 
ual SEDs from \ym to 200 keV (see Section [3~TT ). By neglecting 
the IR bump we avoid counting twice the UV emission repro- 
cessed by dust. The hard X-ray bolometric correction, fcbd, is 
computed as the ratio between Lboi.i^m and the X-ray luminos- 

ity, L[2-10]keV- 

A remarkably good correlation is found between &b i and 
the Eddington ratio. Since the choice of the independent or de- 
pendent variable is not straightforward, we here computed the 
OLS bisector for the fcboi - ^Edd relation as already done for the 
L2500 A ~~ L2 keV relation. We found that the best-fit parameters 
for the fcboi - /lEdd relation using OLS(Y|X) (i.e. treating /l E dd as 
the independent variable) are 

Logk b oi(/l E dd) = (0.273±0.045)Log^ Ed d+(1.656+0.056), (12) 

while the relation using OLS(X|Y) is 

Logkboi(^Edd) = (1.200±0.170)Logi Ed d+(2.600+0.180). (13) 

We then compute the bisector of the two regression lines and we 
find 

Logk b oi(/l E dd) = (0.643±0.043)Log^ Ed d+(2.032+0.053). (14) 

The absence of a correlation is excluded at the 15<x level. 

Figure [TT] shows k^ Q \ as a function of /l E dd for our sub- 
sample of Type 1 AGN and for the sample of local Seyferts by 



Fig. 9. Plot of a ox vs. redshift for our X-ray selected sample: 322 
spectroscopic sources (circles) and 223 sources with photomet- 
ric redshift (crosses). The solid line represents our best-fit given 
by equation (1101) . For comparison, we show also the best-fit de- 
rived by S06 (short-dashed line). The lower panel shows resid- 
uals (Aa ox ) from our best-fit relation, while error bars represent 
the mean and the lcr dispersion of the residuals for each Az = 1 
bin. 



VF09, where radio-loud objects and low X-ray flux observations 
are removed. Points with error bars are obtained averaging the 
values of bolometric correction in each bin of /l E dd, and stan- 
dard errors at lcr are plotted for comparison (see Section 4.3 of 
VF09 for details). The VF09 sample contains AGN with simul- 
taneous optical, UV and X-ray data retrieved from the XMM- 
Newton EPIC-pn and Optical Monitor (OM) archives, while the 
virial Mbh is estimat ed using the reverberation mapping method 
dPeterson et al.ll2004l) . It is worth noting that, although the meth- 
ods to construct the SEDs and measure Mbh are completely dif- 
ferent from those adopted by VF09, the trend of increasing bolo- 
metric correction with Eddington ratio is confirmed, with mean 
(fcboi) ~ 22 for A Edd < 0.1, (fcboi) ~ 27 for 0.1 < A EAd < 0.2 and 
(k boi ) ~ 53 for A Edd > 0.2. 

Equation ( fT2l is in agreement, within the errors, with the 
VF09 results, which are shown in Figure [TT] The equation (fl2l . 
represented by the orange dashed line, and red bins are directly 
comparable, because of the treatment of the independent vari- 
able. 

Differently from VF09, we found a correlation also between 
a ox and the Eddington ratio. In the same interval of A Edd we have 
a larger number of sources, hence the presence of the correlation 
could be simply due to the better statistic, although the disper- 
sion is rather large. In Fig.[T2]we present a m against A Edd , where 
the best-fit relation using OLS(Y|X) is 

ttox^Bdd) = (0.133 ± 0.023)Log A Edd + (1.529 + 0.028), (15) 

while the best-fit relation using OLS(X|Y) is 

a m (A Edd ) = (0.719 ± 0.127)Log A Edd + (2.124 + 0.132). (16) 
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Fig. 10. Plot of the hard X-ray bolometric corrections against 
a ox for the 343 Type 1 AGN with detection in both soft and 
hard bands. The blue solid line represents the best-fit relation 
from eq.dTTb. while the long and the short dashed lines repre- 
sent the lcr and the 3<x dispersion in the distribution of bolomet- 
ric correction, respectively. The bottom panel shows residuals 
(A[Log kboi, l^m]) from the best-fit relation, while error bars rep- 
resent the mean and the lcr dispersion of the residuals for each 
A(Log k b oi, l/im) = 0.2 bin. 



Finally, the bisector of the two regression lines is 

ffox^Edd) = (0.397 + 0.043)Log A Edd + (1.797 ± 0.047). (17) 
The absence of a correlation is excluded at ~ 9<x level. 



5. Effects of reddening and host-galaxy light 

As mentioned in Section [3~T| and discussed at length in Elvis et 
al. (2009, in prep.), the objects used in this analysis show a large 
variety of SEDs. In addiction to objects with a "typical" Type 1 
AGN SED (see Fig. |4]i, there is also a not negligible number of 
objects which show a red optical-UV SED (see Fig. |5). While 
some of them may be intrinsically red AGN, it is likely that for 
most of them this red optical-UV SED is due either to intrinsic 
absorption or to a signif icant contribution o f emission from the 
host-galaxy, or both (see Richards et al. 2003, hereafter R03). In 
the analysis presented in the previous Sections, we used the "ob- 
served" L 2500 ^. If intrinsic absorption or contribution from the 
host-galaxy are not negligible, the used L 2500 ^ would be biased. 
In particular, the intrinsic L 2500 ^ would be higher in presence 
of reddening, while it would be smaller in presence of a sig- 
nificant contribution from the host-galaxy. In the following we 
try to estimate the size of the effect of extinction and/or host- 
galaxy contribution to the relations between L 2500 ^ -L2 kev and 
a ox - L 2500 ^ discussed in Section |4~Tl and |4~2l 



Fig. 11. Hard X-ray bolometric correction versus Eddington ra- 
tio for the 150 Type 1 AGN with BH mass estimate. The solid 
black line shows the best-fit relation that we found using the OLS 
bisector algorithm (see eq. lfT4l ). while the orange and the green 
dashed lines represent equations ( fT2] i and ( fT3l . respectively. Red 
open circles represent the sample by VF09 (25 sources, see their 
Fig. 6 for details) and the corresponding bins with error bars. 




0.001 0.01 0.1 1 



Fig. 12. Plot of a QX as a function of Eddington ratio for the 150 
Type 1 AGN and the corresponding bins with error bars. Key as 
in Figure [TT] The solid black line shows the best-fit relation that 
we found using the OLS bisector algorithm (see eq. 1(1711 ). while 
the orange and the green dashed lines represent equations ( fT"5l ) 
and ( fTol l, respectively. 
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*te*- r > logW^sooA] [^g s- Hz-] 



Fig. 13. Distribution of redshifts versus the relative A(g + - f ) 
color for the sample with z < 4. The dashed lines represent, 
from left to right, the expected loci for E(B-V)=0.04, 0. 12, 0.20, 
respectively. Open circles and squares mark reddened AGN: 
blue corresponds to E(B-V)~0.12 (16 AGN), magenta to E(B- 
V)~0.16 (16 AGN), cyan to E(B-V)~0.20 (14 AGN) and yellow 
to E(B-V)~0.24 (7 AGN). 



5.1. Intrinsic Extinction 

Following the procedure described in R03, we estimate the un- 
derlying continuum color by subtracting the median colors of 
AGN in redshift bins of Az = 0.2 from the observed color of 
each source in that bin. W e refer to A(g + -;'*) - (g + -i*)-(g + -i*) 
as the relative color (see Richards et al. 2001), where the aver- 
age (g + - i*) is computed in each redshift bin. All sources have 
been detected in the Subaru g + band, only 8 sources do not have 
i* CFHT magnitude. For these AGN we consider the detection 
in the i + Subaru band. In Figure [13] we plot the redshift of the 
sources as a function of their relative color (see Fig. 6 in R03). 
The A(g + - i*) distribution shows a large scatter with respect 
to the SDSS sample in R03. Given the quality of the COSMOS 
photometric data, this is not due to photometric errors, but likely 
to the different selection criteria: SDSS quasar candidates are 
selected using optical color-color selection, so the scatter in 
A(g + -i*) is smaller than in our X-ray selected sample. However, 
we can still use the same plot in order to identify possibly red- 
dened sources. Under the assumption that all Type 1 AGN have 
the same continuum shape, the dashed lines show the expected 
change in rela tive color as a fun ction of redshift for an SMC 
reddening law dPrevot et al.lll984h with E(B-V) = 0.04, 0. 12 and 
0.20 moving from left to right in the A(g + - f ) axis. Following 
R03, we define dust-reddened Type 1 AGN all the sources that 
lie to the right of the dashed line at E(B-V)=0.12. Using this 
definition, about 10% of the total sample is affected by intrinsic 
absorption. For about 80% of the sample, reddening is negligi- 
ble, with E(B-V) less than 0.04. The 53 reddened Type 1 AGN 
have been divided in 4 subsamples which correspond to an av- 
erage E(B-V) value of about 0.12, 0.16, 0.20 and 0.24 (see cap- 



Fig. 14. Distribution of L 2500 ^ for the unreddened sample (open 
histogram) and for the 53 sources with E(B-V)> 0.12, before 
(red hatched histogram) and after reddening correction (blue 
filled histogram). 



tion in Fig. [T3Tl. Monochromatic luminosities at 2500 A of these 
objects were corrected using the SMC reddening law and the 
corresponding average value of E(B-V) in each bin. Figure [T4l 
shows the distribution of the optical luminosities before and af- 
ter de-reddening. The L 2500 ^ distribution of the reddened AGN 
is significantly different (i.e. lower luminosity) from that of the 
total sample. This would suggest t hat extinction is mo re impor- 
tant for lower luminosity AGN (see lGaskell e t al. 2004). The av- 
erage shift induced by the correction for the intrinsic extinction 
in the 10% of the total sample is (ALog ]^250(>k 

) = 0.28 + 0.07. 

X-ray absorption is generally negligible in Type 1 AGN; 
however, it is known that a fraction of the order of 10% of broad- 
line AGN may be obscur ed by column densities up to 10 22 cm" 2 
(see lMainieri et al.ll2007[) . Unabsorbed X-ray fluxes can be com- 
puted if the absorption column density (Nh) is known, which 
is not the case for most of the sources in our sample. Hardness 
ratios may be used instead, but they are almost insensitive to col- 
umn densities of the order of 10 22 cirT 2 or slightly higher at the 
average redshifts of the X MM-COSMOS so urces, and they tend 
to over-estimate Nh (e.g. iPerola et"aTll2004l) . In order to quan- 
tify the average impact of X-ray absorption on the a ox distribu- 
tion and bolometric corrections, we have assumed that 10% of 
the sources in our sample are obscured by a column density of 
10 22 crrT 2 . We note that this assumption is likely t o overestimate 
X-ray absorption in Type 1 AGN (see Fig. 13 in lMainieri et al.l 
2007). By correcting X-ray monocromatic fluxes at 2 keV for 
randomly chosen 10% of the sources, the unobscured X-ray 
fluxes are 10% higher. The shift induced by this correction in 
the 10% of the total sample is (ALog L 2 k e v) = 0.04. 

Broad absorption- line quasars (BAL QSOs) are known to be 
X-ray obscured (e .g., iGreen et al.1 1 19951 (Gallagher et al.lll999l 
iBrandt et al.ll2.000l) . and are not included in previous studies of 
optically selected samples because they can cause an artificial 
steepening of the a ox - L 2500 ^ and L 2500 ^ - L2 k e v correlations. 
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Fig. 15. Distribution of the optical index for the total sample 
(open histogram), for the 145 AGN which have a host-galaxy 
contribution less than 10% (filled histogram) and for the 96 
AGN which have a host-galaxy contribution greater than 10% 
(hatched histogram). 



Due to the lack of a systematic analysis of the optical spec- 
tra of broad-line AGN in the zCOSMOS and IMACS spectro- 
scopic surveys and the inclusion of photometric Type 1 AGN, 
we do not have an estimate of the BAL fraction in our sam- 
ple. Assuming they c onstitute about 10-15% of the entire quasar 
population (e.g. R03. lHewett & Foldl20"03l) . BAL QSO may be 
numerous among red sources in Fig[14] and/or X-ray obscured 
AGN. Therefore, we expect that the effects of dust reddening and 
X-ray absorption considered and quantified in the previous para- 
graphs take into account, at least in a statistical sense, the BAL 
QSO contamination. Moreover, the considered fraction of BAL 
QSOs comes from optically selected samples, so that the same 
fraction should be smaller (and possibly redshift dependent) for 
soft X-ray selected samples. 

Applying these corrections for absorption in the optical and 
in the X-ray, we find that the slopes of the L 2500 A - L2 kev and 
a ox -L 25Q0 x relations become /3 corr = 0.782 + 0.021 (OLS bisec- 
tor algorithm) and 0.147 + 0.008. Both these "corrected" slopes 
are within one sigma from those derived with no correction for 
absorption (see eq. J6) and eq. J7)). The slight steepening of the 
^2500 A _ L2 kev relation is due to the fact that the fraction of op- 
tically extincted AGN is higher at lower optical luminosity (see 
Fig. fl~4T>. The corrections for absorption increase the mean value 
of a ox by only 0.01. We therefore conclude that the absorption 
corrections do not change significantly our overall results. 

5.2. Host-Galaxy Contamination 

In order to estimate the possible effect of the contribution of 
the host-galaxy stellar light to the 2500 A luminosity, we have 
considered a subsample of 241 sources for which we have a 
host-galaxy contribution computed subtracting a point-spread- 
function scaled to a central 4 pixel aperture, and measuring the 



exten ded flux in ACS F814W dJahnke et al.ll2004l Uahnke et al.1 
120091) : 145 sources have an extended host-galaxy contribution 
of less than 10%, while for 96 AGN the contamination of the 
host-galaxy light is greater than 10%. This is a robust, model 
independent measure to determine whether substantial extended 
flux is present or not. A useful parameter that we can compute 
for the total sample is the optical spectral index, o-opt. By linking 
it with the morphological information we can try to quantify the 
host-galaxy contribution for each AGN. 

Following lVanden Berk et all (1200 ll) . we compute the optical 
spectral index using the monochromatic luminosities at 2000A 
and 5000A 

a opt = -2.5 Log . (18) 

L 5000A 

Figure [15] shows the distribution of ar opt for the total sample, for 
the 145 sources for which the stellar light contribution to the 
total flux is less than ~ 10% and for the 96 sources for which 
this contribution is greater than 10%. These two samples have 
significantly different distributions of ar opt . Unresolved quasars 
have an average ff p t of ~ 0.5, which is consistent with previous 
results from the SDSS samples, while the sample with a signifi- 
cant galaxy contribution has a range of aopt from to ~ 3 (typical 
of a spiral Sc template) with an average of ~ 1.2. 

We then tried to estimate the fraction of galaxy light at 
2500 A assuming a typical slope a opt = 0.5 for the nucleus com- 
ponent and ff opt = 3 for the galaxies. We further assume that the 
monochromatic luminosity in the range 2000 - 5000A is due to 
the contribution of two power laws: 

L v = Av-° 5 +Gv~\ (19) 

where A and G are the normalization constants we want to esti- 
mate. With this assumption the observed a opt is a function only 
of the ratio between A and G, therefore, at any given a opt value, 
the expected relative contribution of the nuclear and host-galaxy 
emission can be computed. We divided the total sample in 5 bins 
of aopt, with Affopt ~ 0.5, and corrected the observed L 2500A 
using an average value of A/G in each bin. The shift in the 
monochromatic luminosities at 2500 A induced using this cor- 
rection is (ALog L 250()A ) = -0.15 + 0.01, with a dispersion of 
0.21. 

Applying this correction, we find that the slopes of the 
^2500 A _ L2kev and a ox - L 2500A become 0.660 + 0.022 (OLS 
bisector algorithm) and 0.197 + 0.011, respectively. Both these 
slopes corrected for the possible host-galaxy contribution are 
different at about 3<x from those derived with no correction (see 
eq. |6) and 17)). In this case, the flattening of the L 25()0 A — L2 kev 
relation is due to the fact that ffopt tends to be steeper (i.e. with a 
possibly higher contribution from the host-galaxy) for lower lu- 
minosity AGN. This correction decreases the mean value of ff ox 
by ~ 0.06. 

At face value, the possible effects of host-galaxy contribu- 
tion estimated above are not negligible. However, we stress that 
these estimates have to be taken as upper limits, for at least two 
reasons. First, in our derivation of the constants A and G we 
have assumed that all AGN have the same intrinsic a opt ; in pres- 
ence of a dispersion of a opt distribution of AGN, the contribution 
from the host-galaxy would be smaller than that we have derived. 
Secondly, the same would happen also if, as it is likely to be the 
case, both effects (extinction and host-galaxy contribution) are 
at work at the same time. 
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6. Discussion 

6.1. Sample biases and systematics 

We used the COSMOS multi-wavelength database to build a 
large X-ray selected sample of Type 1 AGN and study their op- 
tical to X-ray properties. The sample size is comparable to that 
of op tically selected sam ples previously reported in the literature 
(e.g., Gibson et al. 2008), but to our knowledge this is the largest 
complete X-ray selected sample for which the study of the a ox 
distribution is performed. In order to keep the selection criteria 
as clean as possible and to cope with the lack of spectroscopic 
information at faint optical magnitudes, the Type 1 classification 
is based on both spectroscopic and photometric redshifts. This 
choice guarantees a well defined and relatively simple selection 
function, but it is not completely free from other biases. In par- 
ticular, contamination from photometrically misclassified Type 
1 AGN cannot be excluded (see Section l2T2l >. Other possible bi- 
ases and systematics include the contribution of the effect of dust 
reddening and gas absorption at UV and X-ray frequencies and 
the host-galaxy light to the 2500 A nuclear flux. Even though we 
tried to estimate and quantify the impact of the above mentioned 
biases (see Section 15.11 and I5.2t . some residual contamination 
due to either one or more of these biases cannot be completely 
excluded. However, we are confident that the most important re- 
sults of our analysis, and especially those concerning the average 
properties of the total sample, are not significantly affected. 

The correlations between optical and X-ray luminosities ob- 
tained from the analysis of the spectroscopic sample are recov- 
ered, if only photo-z are considered. Even though the best fit pa- 
rameters of the correlations are slightly different, the inclusion 
of photo-z sources allow us to extend the study of optical and 
X-ray properties to much lower luminosities. 

6.2. Constraints on the X-ray emission models 

The optical to X-ray spectral index a ox connects the two por- 
tions of the AGN broad-band spectrum dominated by the accre- 
tion power and thus it is expected to be a reliable tracer of the 
accretion properties and, in particular, of the relation between 
disk emission, peaking in the UV, and coronal X-ray emission. 
Previous studies have shown a non linear dependence of the 
X-ray versus UV luminosity, which implies a correlation be- 
tween the optical-UV to X-ray luminosity ratio on monochro- 
matic luminosity and/or redshift. Based on exte nsive analysis 
of large samples of optic ally selected AGN (i.e. IVignali et al.l 
I2003L IStrateva et al.ll200l S06, J07) it is concluded that a ox is 
primary dependent on optical luminosity at 2500 A. The results 
presented in this paper confirm and extend the previous findings 
to a large sample of X-ray selected AGN, suggesting that band 
selection does not significantly modify the current observational 
picture. The observed L 2500 ^ - L2 kev (or a ox - L 250 o a) corre- 
lation implies that more optical luminous AGN emit less X-rays 
per unit UV luminosity than less luminous AGN. The predicted 
UV to X-ray luminosity ratio depends on the hot corona cover- 
ing factor, optical depth and electron temperatures. The optical- 
UV and X-ray emission of radio-quiet AGN can be explained 
by the interplay between hot electrons in a coronal gas and a 
colder accretion flow. Soft photons from the accretion disk are 
Comptonized by hot electrons and lead to the formation of a 
power law spectrum in the hard X-rays acc ompanied by a high 
energy cut -off at the electrons' temperature (Ha ardt & Maraschil 
Il991ll l993h . If only a fraction of the accretion pow er is released 
in the hot phase, as in the patchy corona model (Haard fet alJ 



1994), the L\jy/Lx ratio is higher than the value computed us- 
ing a model with more uniform corona. The observed correlation 
suggest that disk-corona parameters are depending on UV lumi- 
nosity. 

6.3. Using a ox as a bolometric correction estimator 

An accurate determination of bolometric quantities is essential 
to compute AGN accretion rates and, more in general, to all the 
arguments related to the accretion onto SMBH at large. These 
parameters require to sample the AGN SED over a broad range 
of wavelengths from near-infrared to hard X-rays. Owing to the 
excellent quality of the multi-wavelength COSMOS database, 
this is becoming possible for large AGN samples (Elvis et al. 
2009). Because the optical-UV to X-ray portion of the spectrum 
contains about 60% of the total bolometric luminosity, we have 
investigated to what extent the UV to X-ray luminosity ratio can 
be considered a reliable proxy of the bolometric correction. 

We computed bolometric corrections from the 2-10 keV en- 
ergy range for all the objects in our sample with detection in both 
soft and hard bands. Given the strong luminosity dependence of 
a ox upon UV luminosity which is close to the SED maximum, 
it is not surprising to find a significant correlation between a ox 
and fcboi- However, it is important to note that the correlation, 
best-fitted by a second order polynomial, has a small dispersion 
around the best-fit (see Fig. [TOt , The existence of such a tight 
relation between these two parameters suggests that the AGN 
bolometric output is well traced by a ox over a broad range of 
redshifts and luminosities. Moreover, it offers the opportunity 
to estimate reliable bolometric corrections, for large AGN sam- 
ples, provided that rest-frame 2500A and 2 keV luminosities are 
known. 

6.4. Bolometric corrections and Eddington ratios 

Accretion rates are then estimated using the previously com- 
puted bolometric corrections. A correlation is found between the 
2-10 keV bolometric correction and the Eddington ratio for a 
sizable subsample of 150 objects (see Fig. ITTh. for which black 
hole masses are computed from the broad em ission lines FWHM 
(see Merloni et al. 2009. iTrump et alj|2009l for a detailed dis- 
cussion). The presence of this correl ation was originally sug- 
gested by IVasudevan & FabiarJ d2007l) using broad-band obser- 
vations of nearby Seyfert galaxies and confirmed in a subse- 
quent paper ( VF09) using simultan eous optical-UV and X-ray 
data. A similar trend is reported by Bianc hiet al.l (12009) from 
a multi-wavelength analysis of archival XMM-Newton observa- 
tions. Although we do not have simultaneous data or Mbh es- 
timate from reverberation mapping, our findings are consistent 
with the &boi - J-Edd relation by VF09. This correlation may, in 
principle, be induced by the fact that both parameters depend 
on the bolometric luminosity. However, if we look for the effect 
of the bolometric luminosity using the Kendall-r partial corre- 
lation analysis, we found that the correlation is still significant. 
This suggests that the systematic effects introduced by the pres- 
ence of bolometric luminosity on both axes are not relevant (see 
also Section 5.1.4 in VF07). 

A fairly significant correlation, albeit with a large scatter, is 
also found between a ox and the Eddington ratio. In the frame- 
work of the disk-corona models discussed above, sources with 
steep a ox are explained by a patchy corona where the number of 
blobs or their covering factor decreases for increasing values of 
a ox . For a given BH mass, bright UV emission is due to a high 
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accretion rate. If in highly accreting sources most of the power 
is dissipated in the disk, rather than in the hot corona, then a 
correlation between a ox and Eddington ratio is expected. 

7. Summary and conclusions 

In this paper we have presented the analysis of the UV to X-ray 
properties of 545 radio-quiet X-ray selected Type 1 AGN, in the 
multi-wavelength COSMOS survey. The full data-set covers a 
large range of redshifts (0.04 < z < 4.25) and X-ray luminosi- 
ties (40.6 < Log Lp-io]keV ^ 45.3). It is mostly composed by 
spectroscopically selected Type 1 AGN (322 sources, ~ 60% of 
the total sample), with the addition of 223 Type 1 AGN classi- 
fied on the basis of the best-fitting SED procedure provided by 
S09. We constructed single SEDs for the full sample and, from 
these, we computed optical and X-ray rest-frame luminosities 
at 2500 A and 2 keV, respectively; we also investigated the de- 
pendence of a ox upon redshift, L2 k e v and L 2500 \ using the fully 
parametric EM (estimated and maximized) regression algorithm 
and Kendall-r partial correlation analysis. Moreover, we esti- 
mated bolometric correction and bolometric luminosities for a 
subsample of 343 AGN for which we have detections in both 
soft and hard bands. Our principal results are the following: 

1 . We confirm the L2 kev - L 2500 a correlation parametrized by 
L 2 kev oc if „ , where B = 0.760 + 0.022. The best-fit slope 

ZKev 2500 A' h f 

agrees with previous studies based on optically selected sam- 
ples, which found a value of /3 inconsistent with unity. 

2. The mean value of a ox for the full X-ray selected sample is 
(a ox ) ~ 1.37 + 0.01 with a dispersion around the mean of 
0.18. We confirm the correlation between a ox and L 2500 ^ at 
the 17 cr significance level. The correlation becomes stronger 
if we take into account the effect of redshift using the partial 
correlation analysis (~ 21 cr). The slope of our the best-fit 
relation between a ox and L 2500 ^ is consistent at least within 
~ 1 .6cr with that obtained for an optically selected sample. 

3. From the EM regression method, we find a weak correlation 
between ff ox and redshift; if the effect of L 2500 - A is taken into 
account, the correlation disappears. 

4. We do not find any significant correlation between ff ox and 
L2kev (~ 1-5 cr significance level if we take into account the 
effect of the optical luminosity). 

5. We found a tight correlation between the a ox values and the 
hard X-ray bolometric correction values. The ff ox - k\, a i rela- 
tion could be used as a practical tool to provide an accurate 
estimate (~ 20% at lcr) of the bolometric correction using 
only the a ox value. 

6. We found a correlation for both ff ox and feboi with Eddington 
ratio. Our results for the feboi - ^Edd relation and for a ox - ^dd 
suggest that there is a connection between the broad-band 
emission, mostly in the UV, and the Eddington ratio, which 
is directly linked to the ratio between mass accretion rate, 
M acc , and Eddington accretion rate, M Edd . 

7. We have estimated the possible effects of intrinsic absorption 
and/or contribution of the host-galaxy on the L 2500 ^ - L2 kev 
and a ox - L 2500 ^ relations. If the main reason for the red 
optical-UV SEDs in our sample were intrinsic extinction, the 
resulting L 2500 ^ - L2 kev relation would be slightly steeper 
(B = 0.782 + 0.021); if, instead, the host-galaxy contri- 
bution were important in determining the red optical-UV 
SEDs, the resulting L 2500 ^ - L 2 kev relation would be flat- 
ter (J3 = 0.660 + 0.022). Neither effect can produce a linear 
correlation (j3 = 1). 
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